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^ , Abstract 

Motivated by recent observations which detect dust at large galactocentric distances 

I in the disks of spiral galaxies, we propose a mechanism of outward radial transport 

' of dust by spiral stellar density waves. We consider spiral galaxies in which most of 

■ dust formation is localized inside the corotation radius. We show that in the disks 
I of such spiral galaxies, the dust grains can travel over radial distances that exceed 
' the corotation radius by roughly 25%. A fraction of the dust grains can be trapped 
. on kidney-shaped stable orbits between the stellar spiral arms and thus can escape 

^ i the destructive effect of supernova explosions. These grains form diffuse dusty spiral 

I arms, which stretch 4-5 kpc from the sites of active star formation. About 10% of 

I ' dust by mass injected inside corotation, can be transported over radial distances 

2 , 3 — 4 kpc during ~ 1.0 Gyr. This is roughly an order of magnitude more efficient 

' than can be provided by the turbulent motions. 
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1 Introduction 



Interstellar dust has been recently recognized to spread over large distances 
from their source, and thus can be considered as a tracer of dynamical pro- 
cesses responsible for the circulation of material in disk galaxies. Observa- 
tions of edge-on galaxies reveal a great amount of dust far outside the galactic 
planes: in NGC 891 and NGC 4212 dusty clouds extend up to 2; ~ 2 kpc above 
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the plane, and contain in total as much as ~ 10^ Mq (Howk and Savage, 1997, 
1999; Rossa et al., 2004, Thompson et al., 2004). Several mechanisms can con- 
tribute to transport of dust in the vertical direction: convective motions or 
bore flows (Gomez and Cox, 2002) associated with spiral density waves, chim- 
neys produced by multiple supernova explosions (Norman and Ikeuchi, 1989), 
and radiation pressure (Ferrara et al., 1991; Shustov and Vibe, 1995; Dettmar 
et al., 2005). 



The situation is less clear when the radial distribution of dust is considered. 
Observations of spiral galaxies in the far infrared band show that dust can ex- 
tend in the radial direction far beyond the stellar disks. Neininger et al. (1996) 
found that the radial distribution of 1.2 mm dust emission from NGC 4565 is 
two times wider than CO emission. ISO observations of several spiral galaxies 
at 200 /im (Trewhclla ct al., 2000; Bianchi ct al., 2000) revealed the presence 
of relatively cold dust outside the stellar disks with a scale length Rd ^ 1.5 
{Rd and R* are the radial scales of dust and stars, respectively). A similar 
conclusion was made by Xilouris et al. (1999) for four spirals observed in the 
optical and near infrared bands. Bianchi et al. (2000) found that the best fit 
of the 200 //m emission for NGC 6946 gives Rd ^ 3 R^. MHD turbulence can, 
in principle, spread dust outward in radial direction (Cho ct al., 2003), how- 
ever, the corresponding time for galactic scales seems to be too long. Indeed, 
the coefficient of turbulent diffusion D ~ {Lv)/3 is of the order of 10^^ cm^ 
s~^ for typical interstellar values L ~ 100 pc and i> ~ 10 km s~^. Therefore, 
dust particles can diffuse by turbulent motions over scales ^/{AR^ > 3 kpc 
in a characteristic time t > 16 Gyr. In these conditions, only regular hydrody- 
namic motions can explain the presence of dust outside stellar disks in spiral 
galaxies. 



In this paper we explore the possibility that the radial transport of dust in 
spiral galaxies is driven by the hydrodynamic motions associated with spiral 
stellar density waves. In Sect. 2 a set of model hydrodynamic equations for 
dusty interstellar medium is described and the equilibrium parameters are 
fixed. In Sect. 3 the generation of spiral pattern in the stellar disk is formulated. 
Sect. 4 describes trajectories of single dust particles and Sect. 5 describes the 
dynamics of dust fluid as a whole including formation of dusty spiral arms. 
Sect. 6 summarizes our results. 
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2 Model description 



2.1 Basic equations 



Our model disk is a two-fluid system composed of dust and gas, which evolves 
in the external potential of the stellar disk and dark matter halo. We assume 
that the dust and gas are coupled by friction. The gas disk is isothermal at T = 
10^ K and the dust pressure is negligible. We use the thin disk approximation, 
i.e. all quantities are vertically integrated. The system is described by the basic 
hydrodynamic equations, where we denote the dust and gas by the subscripts 
d and g, respectively. 

The continuity equation for the gas and dust is 



dt 



+ V-(Sg,dt;g,d) = 0, (1) 



where Eg_d are the gas and dust surface densities and Vd,g are their velocities. 
The equations of motion are 



dv„ , „^ VPg Ed/ 

-gf + K- V)t;g = - V$,i,,2,h - ^ + (2) 

-^ + ivd- V)t;d = - V$si,s2,h - /, (3) 

where, $si!*^'s2i and $h are the contributions to the total gravitational po- 
tential by the axisymmctric stellar disk, non-axisymmetric spiral density per- 
turbation, and dark matter halo, respectively. The friction force is defined as 
/ ^ A(vd - Vg) (see Draine and Salpeter, 1977; Noh et al., 1991). 

We take the friction coefficient A to be a simple function of the disk conditions. 
We assume that the thickness of the dust disk is comparable to that of the 
gas disk. Then the coUisional timescale Tc of energy and momentum exchange 
between the gas and dust particles may be used as an estimate of A~^, 

A ^ = WWn,^ 
rrid 



where rig is the number density of gas particles, rUg and rud are the masses of 



gas and dust particles, and = ^JSRT/jj, is the gas velocity dispersion; t^. is 
estimated for the subsonic relative motion of neutral dust grains and gas. Since 
we use the thin disk approximation, the gas volume density in Eq. (4) is ap- 
proximated assuming a local vertical hydrostatic equilibrium mgUg = Sg/^o, 
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where the gas disk thickness is defined as^ zq — 2Cg/(7rGEg). where Cg is the 
sound speed. The geometrical cross section is used as an estimate of the coUi- 
sional cross section Uc = tto^, where a is the average radius of a dust particle. 
In the present simulations, we assume rrid — 10~^^ g and a — 10~^ cm. It 
is well known that interstellar dust grains are charged, and Coulomb inter- 
actions with charged ambient ions and electrons increase the friction force 
(Draine and Salpeter, 1979; Weingartner and Draine, 2001). However, as we 
are mostly concerned with the diffuse HI phase of the ISM (warm neutral 
medium, WNM, in standard nomenclature), this increase does not crucially 
change the resulting friction coefficient. Indeed, the Coulomb friction can be 
written as Fq = x(Z)F^j, where F^j is the friction in an HI environment, 
x(Z) = (a^/a)^ In A a;j(mi/mp), az = Ze'^/kT, InA is the Coulomb loga- 

i 

rithm, Xi is the fractional ionization of i-th ions or electrons, rrii is the corre- 
sponding mass, averaging is over charge distribution of dust grains (Weingart- 
ner and Draine, 2001). For typical conditions in the WNM; x ~ 0.1, n ~ 0.3 
cm~^, and dust charge (Z) ~ 60 (Yan et al., 2004) we get x ~ 1. Therefore, 
we will restrict our consideration with the friction coefficient A as determined 
for the neutral dust. 



2.2 Equilibrium gas disk 



Initially, our model galaxy consists of a rotating gas disk balanced by the grav- 
ity of an axially symmetric stellar disk and spherically symmetric dark matter 
(DM) halo. The dust disk is introduced later when the non-axisymmetric part 
of the stellar gravitational potential is set. The rigid stellar disk is assumed to 
have a power-law radial density profile of the form 



Ss(r) 



2^ 



-3/2 



(5) 



where B'^ = 2nGrf'Eso- The gravity force of such a density distribution is given 
by Toomre (1963) 



dr 



-(r^ + r^) 



-3/2 



(6) 



In the following, we use the central stellar density Sgo = 1.2 x 10^ Mq pc ^ 
and Ts = 3 kpc, which gives us a total stellar mass of Mgt = 7 x 10^° Mq. 

^ Since the dominant component of gravity in the vertical direction is the stellar 
disk, the total surface density (Stot = + Sg) should appear in the expression for 
Zq. This may increase the friction coefficient A, especially in the inner regions. We 
discuss this issue in Sect. 5 
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The density distribution of the rigid DM halo is assumed to be that of a 
modified isothermal sphere (Binney and Tremaine, 1987) 

{1 + r/rhy 



where the central density pho and the characteristic scale length are given 
by Mac Low and Ferrara (1999) and Sihch and Tenorio-Tagle (2001): 



— 1/3 

Pho = 6.3 X 10^° ( ^] h-'l' Me kpc"^ (8) 



rh = 0.89 X 10-^ {^\ h}'^ kpc. (9) 



Here, h is the Hubble constant in units of 100 km s^^ Mpc^^ and — 
10^^ Mq is the total DM halo mass. We note that the DM halo mass within 
the computational domain r < 20 kpc is Mh(r < 20 kpc) = 5 x 10^*^ Mq. We 
adopt h = 0.65 throughout the paper. The gravitational force of spherically 
symmetric DM halo can be expressed as 



= 47rG'pho?^h [r/rh - arctan(r/rh)] ( — j . (10) 



The gas disk has an exponentially declining density profile with the central 
surface density Ego = 30 Mq pc~^, and radial scale length — 9 kpc. The 
total mass of the gas disk within the computational domain (r = 20 kpc) is 
Mg = 1.0 X 10"'^'' Mq. The gas disk contains only a small fraction (~ 10%) of 
the total mass in the computational domain. Hence to a first approximation 
we can assume that the gas moves in the potential field of the stars and DM 
halo. 

Once the density profile of the gas disk is fixed, its initial rotation curve 
(RC) is obtained by solving Eq. (2), where the friction force / is set to zero 
and the time variations are neglected. The resulting initial RC is plotted in 
Fig. 1, where it reaches a maximum circular velocity of 200 km s~^ at r ~ 
3 kpc. Beyond this maximum, the rotational velocity undergoes a gradual 
decline. Because the gas surface density has a slowly declining radial profile 
(rg = 9 kpc) and the gas is assumed to be isothermal, the contribution of the 
pressure gradient in the gas rotation curve is less than 1% everywhere in the 
disk. 
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Fig. 1. The initial rotation curve of the gas disk. 
2.3 Computational techniques 

A set of hydrodynamic equations (l)-(3) in polar coordinates (r, 0) is solved 
using the method of finite-differencing with a time-explicit, operator-split so- 
lution procedure (ZEUS-2D) described in detail in Stone and Norman (1992). 
We use a resolution of 200 x 200 grid points in a polar computational domain 
with a radius of 20 kpc. We set a free boundary condition at the outer bound- 
ary, i.e. the gas and dust are allowed to flow out freely from the computational 
domain. The outflowing gas/dust is assumed to be lost by the parent galaxy. 



3 Generation of spiral pattern 



Once the axisymmctric equilibrium gas disk is constructed, wc introduce the 
non-axisymmetric spiral gravitational potential of the stellar disk $s2- The 
latter is deflned in the form of a running spiral density wave (Lin et al., 1969) 

^s2 = -C(r) cos [m(cot(i) ln{r/r^p) + - Q^pt)] . (11) 
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Fig. 2. Behaviour of n±K/2, where k is the epicychc frequency and O is the angular 
velocity of gas. The angular velocity of spiral pattern Ogp is shown by the solid line. 

Here, C(r) is the amplitude of the spiral gravitational potential, m is the 
number of spiral arms, i is the pitch angle, rgp is the characteristic winding 
length of spiral arms, and fisp is the angular velocity of spiral pattern. In the 
following we adopt m = 2, i = 25°, Vgp = 6 kpc, Qsp = 24 km s^^ kpc~^. 
This choice of flgp places the corotation of the gas disk at ~ 8.0 kpc and 
the outer Lindblad resonance at ^ 12.0 kpc. The inner Lindblad resonance 
is absent. The behavior of gas angular velocity fl and Q ± k,/2 is shown in 
Fig. 2. According to Kennicutt (1981), a pitch angle of i = 25° is typical for 
Sc galaxies. The efficiency of radial transport of dust for different pitch angles 
is discussed in Sect. 5. 

The logarithmic spiral defined by Eq. (11) is quickly winding up as r — > 0. 
This may create large gradients of gravitational potential near the origin and 
lead to serious numerical difficulties. For instance, if the amplitude C(r) is 
independent of r, than the ratio of maximum non-axisymmetric perturbing 
force to the total axisymmetric gravity force /3(r) = max\'V^s2\/{d^si/dr + 
d^Yi/dr) sharply increases as r — 0. This is illustrated in Fig. 3, where /?(r) 
is plotted by the dotted line for C = 0.005 (in dimensionless units). The 
apparent growth of f3 with radius for r > 3 kpc is also unrealistic since the 
spiral arms are expected to vanish at larger radii. In view of these difficulties, 
we use a radially dependent amplitude of the form C(r) = Co{r)°'^'^\ Here, 
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Fig. 3. The ratio of the maximum, non-axisymmetric perturbing force 
|V<I>s2| = [(9$s2/9r)^ + {r^^d^s2/dcl))'^Y^'^ to the total axisymmetric gravity force 
(d^si/dr + d^i^/dr) as a function of radial distance. The solid and dashed lines are 
explained in greater detail in the text. 

Co(r) is a linear function of r which has a value of at r = kpc (ensuring 
that |V$s2| diminishes as r — >■ 0) and attains its maximum value of 0.0045 
(in dimensionless units) at r = 20 kpc. The exponent a{r) decreases linearly 
with radius from a = 2 at r=0 kpc to a = —2.1 at r=20 kpc. The resulting 
profile of /5(r) is shown in Fig. 3 by the solid line. Since the non-axisymmetric 
gravity force scales as |V$s2| oc l/r^^"*^*'\ the ratio (3 increases with radius 
at r < 4 kpc, and decreases at r ;^ 4 kpc. At the position of corotation 
Tcr ~ 8 kpc, (3 drops by roughly a factor of 2 as compared to its maximum 
value at r f=:i 4 kpc and continues to decrease at larger radii. Wc note that the 
maximum, non-axisymmetric perturbing force never exceeds 12% of the total 
axisymmetric gravity force. 

The non-axisymmetric part of the stellar gravitational potential $s2 is turned 
on slowly. Specifically, $s2 is multiplied by a function e{t), which has a value 

of at t=0 and linearly grows to its maximum value of 1.0 at t > 200 Myr. 
It takes a few hundred Myr for the gas disk to adjust to the spiral distortion 
and develop a spiral structure as shown in Fig. 4. The gas spiral pattern (as 
well as the stellar one) rotates counterclockwise. The position of corotation 
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Fig. 4. Gas density distribution at t = 0.8 Gyr from the beginning of numerical sim- 
ulations. The corotation radius is shown by the circle. The scale bar is in Mq pc~^. 
The contour line sketches the region of supercritical gas density (see Sect. 5 for 
details) . 

is shown by the circle. The strongest gas response to the underlying stellar 
spiral density wave is seen inside the corotation circle. However, a weak spiral 
structure can still be traced outside corotation at radii r = 10 — 15 kpc. 



4 Trajectories of dust particles 

In this section we study the trajectories of single dust particles in the com- 
bined gravitational potential of stellar disk and dark matter halo including 
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Fig. 5. Polar plot (r, 9) showing trajectories of three test dust particles. The radial 
distance is in kiloparsecs and the azimuthal angle is in radians. Stellar spiral arms 
are located along azimuthal angles of and tt. The direction of motion of test 
particles in the frame of reference of stellar spiral arms is indicated by arrows. 

the effects of dust-gas friction. The equation of motion of the dust fluid (3) 
can be obviously reduced to a set of equations of motion for a single dust par- 
ticle and subsequently solved using the Runge-Kutta scheme. We introduce 
test dust particles at different radial distances from the galactic centre and 
at different azimuthal angles 6^ with respect to the position of stellar spiral 
arms. In our simulations, the stellar spiral arms are modeled by the spiral 
gravitational potential rotating counterclockwise in the laboratory frame of 
reference. Therefore, the position of stellar spiral arms at any given time is 
determined as the position of minimum in the spiral gravitational potential 
defined by Eq. (11). All dust particles are injected at t = 0.8 Gyr when the gas 
spiral structure is fully developed (see Fig. 4), and they are assigned velocities 
drawn from the parent gas. We followed the trajectories of test dust particles 
for 2 Gyr. 

The phase (ra, 9d) trajectories of three representative dust particles in the 
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Fig. 6. Schematic description of gravity forces (shown by the dashed arrows) acting 
on test dust particles from the stellar spiral arms (plotted by thick solid and dashed 
lines) . Velocities of dust particles in the frame of reference of stellar spiral arms are 
indicated by solid arrows. Two test particles inside corotation and four test particles 
outside corotation are shown. Arrows are not to scale. 

frame of reference of stellar spiral arms are shown in Fig. 5. The stellar spiral 
arms in Fig. 5 are always located along azimuthal angles ^sp — and tt. 
The dashed hne shows the trajectory of a dust particle placed initially at 
rd=5.0 kpc in the arm (^d = 0; hereafter, particle 1), whereas the solid line 
gives the trajectory of a dust particle placed initially at rd=7.0 kpc in the 
interarm region [9^ = 7r/4; hereafter, particle 2). Dust particles similar to 
particle 1 may be injected by supernova explosions which tend to occur near 
spiral arms. Particle 2 may be created by red giant stars. 

It is apparent that the trajectories of particles 1 and 2 are quite different. 
Particle 1 moves on an elliptical orbit in one direction (counterclockwise), re- 
maining always inside the corotation circle (rcr ~ 8.0 kpc). The trajectory 
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of particle 1 is fairly stable. It comes closer to the origin when traveling be- 
tween the arms (^d ~ ^/2) and it passes through the spiral arms 4 times 
during 0.65 Gyr. Since dust particles are thought to be destroyed by super- 
nova explosions, such a frequent encounter with the spiral arms implies that 
particle 1 may not hve longer than 0.5 Gyr (although it depends on the cov- 
ering factor of sufficiently strong shock waves, vsn ^ 100 km s~^, generated 
by supernova explosions (SNe) and propagating perpendicular to the spiral 
arms). Particle 2 has a stable, kidney-shaped trajectory, which shows an im- 
pressive radial migration of ~ 4 kpc, as opposed to only 1 kpc for particle 1. 
Interestingly, particle 2 is caught between the spiral arms, periodically passing 
through the corotation circle from the inner part of the disk to the outer part 
and vice versa. Since particle 2 never encounters the spiral arms during 2 Gyr, 
it may live significantly longer than particle 1. Dust particles that are born 
at 6'd ~ 7r/2 (i.e. exactly between the spiral arms) but are sufficiently away 
from the galactic center, move on much more compact orbits as shown by the 
dotted line in Fig. 5. 

To summarize, we find that all particles born at r < 6 kpc, i.e. ^ 2.0 kpc 
away from the position of corotation radius Vcr ~ 8.0 kpc, always stay inside 
corotation. This is irrespective of their azimuthal angles of birth. Trajectories 
of those dust particles resemble that of particle 1 shown in Fig. 5 by the dashed 
line. Such particles often encounter spiral arms and thus may be destroyed in 
one rotation period by shocks from SNe. The fraction of surviving dust grains is 
determined by the probability for them to meet a sufficiently strong [vs ^100 
km s~^) destroying shock wave inside the arm. On the other hand, particles 
injected near corotation at 6 kpc <; r < 8.0 kpc (and sufficiently away from 
spiral arms) spend most of their lifetime in the interarm region. Therefore, 
most of these dust particles can survive over the entire migration period in 
the outer galactic regions. 

In the absence of spiral arms, dust particles move on circular orbits deter- 
mined by the axisymmetric global gravitational potential of a galaxy and the 
drag force from the gas disk. When the non- axisymmetric spiral gravitational 
potential is introduced, dust particles move on various non-circular orbits as 
shown in Fig. 5. In order to understand the origin of such non-circular orbits, 
we consider the influence of a non-axisymmetric spiral gravity force on the 
motion of dust particles. The gravity force of spiral arms acting on moving 
particles is schematically plotted in Fig. 6 by the dashed arrows. Six parti- 
cles are presented; two particles inside the corotation circle and four particles 
outside the corotation circle. The direction of motion of dust particles in the 
frame of reference of the spiral pattern is shown by the solid arrows. Two 
spiral arms are schematically shown by the solid and dashed lines. A moving 
dust particle may see the nearest spiral arm as either convex or concave. If 
the dust particle is closer to the concave spiral, the latter acts to increase 
(on average) the radial distance of the particle due to its gravitational drag. 
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Fig. 7. Gray-scale images of the gas distribution (the upper panels), dust distribution 
(the middle panels) and the dust-to-gas ratio (lower panels). The elapsed time is 
indicated in each panel. The position of corotation of the gas disk is outlined by 
the circle. The scale bars in the upper and middle panels give the column density 
in Mq pc~^ 
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On the contrary, if the dust particle is closer to the convex spiral, the latter 
tends to decrease the radial distance of the particle. If we now recall that the 
trailing spiral arms in Fig. 5 are located along the azimuthal angles and 
TT and rotate counterclockwise, then the effect of gravitational drag from the 
spiral arms becomes clearly evident in the trajectories of dust particles. For 
instance, when dust particles arc in the first quadrant, they are attracted by 
the convex spiral located at ^^sp = 0. As a consequence, the radial position of 
these particles decreases. When dust particles enter the second quadrant, their 
motion is controlled by the concave spiral located at = tt and their radial 
position increases. Therefore, the trajectories of dust grains may be sensitive 
to the pitch angle of stellar spiral arms. The maximum radial migrations of 
dust grains are expected for a pitch angle of 45°. 

The action of the friction force is more difficult to assess. The friction between 
dust and gas may help reverse the direction of motion of a dust particle (with 
respect to the spiral arms) when it passes through corotation. Indeed, when 
the dust particle is inside the corotation circle, its velocity is larger than that 
of spiral arms (the same is true for the gas in general). After passing through 
corotation, the dust particle enters the region which is characterized by a 
smaller velocity than that of spiral arms. Therefore, we may expect that the 
friction between gas and dust would act to equalize the velocities of dust 
and gas, and thus decelerate the dust particles. This would eventually reverse 
their direction of motion with respect to the spiral arms as seen in Fig. 5 for 
particle 2. 



5 Response of dust to the stellar spiral density wave 



In this section we study how dust responds to the imposed spiral disturbance. 
The dominant sources of dust formation in disk galaxies arc stellar winds from 
late-type giant and supergiant stars, supernova explosions, and growth of pre- 
existing grains in dense n ^ 10^ cm~^ molecular clouds (Dwek, 1998). All 
of these sources are mainly located inside the corotation circle. Indeed, the 
strongest gas response to the stellar spiral density wave in Fig. 4 is obviously 
seen inside corotation. We compute the critical density for star formation as 
defined by the Toomre criterion (Kennicutt, 1989) 



ipToomre _ '^s 

^crit — _ ^ > V-I-^J 
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where k is the epicychc frequency, and the shear criterion (Martin and Ken- 
nicutt, 2001) 



V shear _ 2^5Ash/i;s 



where Agh = —0.5 r dQ/dr is the local shear rate (the Oort constant). The 
contour line in Fig. 4 traces the region of supercritical density with Eg > 
min(EjJJ;°^°™''®, E^f^'jf'^'') which is evidently localized inside corotation. Although 
the gas density outside corotation may occasionally become supercritical (not 
shown in Fig. 4) especially along the outer diffuse gas spirals, a lower gas 
density and higher temperature as compared to the inner region can work 
against star formation in that region. We thus conclude that most of star 
formation in our model galaxy is localized within the corotation circle. The 
radial scale length of our old stellar disk is rg = 3 kpc which also implies that 
most of low to intermediate mass stars that are capable of producing dust are 
confined within the corotation circle. 

We instantaneously inject 3.4 x 10^ Mq of dust mass inside the corotation 
circle (r < 7.8 kpc) and then determine the fraction of dust that can move 
outside corotation. In the present simulations, the dust particles are injected 
at t = Myr when a spiral structure in the gas disk has not yet developed. 
We also tried to introduce dust at t=800 Myr (when the gas spirals are fully 
developed) and confirmed that the time of dust injection has little influence 
on dust dynamics. The spatial distribution of injected dust is such that the 
dust-to-gas ratio 7 = 0.01 inside corotation, and 7 = outside corotation. The 
velocity of injected dust is equal to that of the gas. A conservative estimate 
of dust-to-gas mass ratio in the solar neighborhood gives an average value of 
7 = 0.006 ranging from 0.002 up to 0.04 in H2 regions (Spitzer, 1978). 

The upper and middle panels in Fig. 7 show two snapshots of the evolution of 
the gas and dust surface densities, respectively, whereas the lower panels give 
the dust-to-gas mass ratio. The numbers in each panel indicate the elapsed 
time since the beginning of the simulation. The circle indicates the position 
of corotation of the gas disk. A strong spiral response to the stellar density 
wave develops in the dust disk after approximately 3 revolution periods of our 
model galaxy. The dust density closely follows that of gas in the inner regions 
of our model galaxy. Small radial and azimuthal variations of the dust-to-gas 
ratio around the initial value 0.01 are indicative of a dominant role of frictional 
force in the dynamics of dust in the inner 5 kpc. This is in agreement with the 
results of Sect. 4 where we have found that dust particles born in the inner 
5 kpc move on fairly stable elliptical orbits around the galactic center. At 
larger radii 6 kpc ^ r < 8.0 kpc, however, significant variations of 7 become 
evident. There are regions near corotation that are completely devoid of dust, 
even though the initial dust-to-gas ratio is 0.01. A noticeable portion of dust is 
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Fig. 8. a) Dust mass ouside corotation as a percentage of the total dust mass. 
The solid and dashed lines correspond to the models with pitch angles of stel- 
lar spiral density wave i = 25° and i = 15°, respectively. The dotted line gives 
the dust mass outside corotation for a lower rotation frequency of stellar spirals 
n^p = 21 km s-i kpc-^ b) G as mass inside corotation as a percentage of the initial 
total gas mass. The solid, dashed, and dotted lines correspond to the same models 
as in the upper panel. The difference in the values of the gas fraction inside corota- 
tion for Qsp 



21 km s ^ kpc ^ and Qsn = 24 km s ^ kpc ^ is due to the different 



positions of the corotation radius in the two cases. 
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observed outside corotation. This dust tends to form into diffuse outer spiral 
arms, which appear to be a smooth extension of the inner dust spiral arms and 
can be traced up to r = 10 kpc. The dust-to-gas ratio in the outer dust arms is 
not dissimilar to that of the inner r < 5 kpc region. We thus conclude that the 
perturbing force of stellar spirals starts to play an important role in the dust 
dynamics near corotation resulting in an outward radial transport of dust. If 
we now identify the regions of supercritical gas density by contour lines in the 
upper panels of Fig. 7, than it becomes evident that the dust component can 
be found as far as 4-5 kpc from the sites of active star formation. 

The solid line in Fig. 8a gives the dust mass outside corotation (as a percentage 
of the total dust mass) versus time elapsed since the injection of dust. Approx- 
imately 10% of the initial injected dust mass is transported outside corotation 
during 1 Gyr. This implies that the dust transport rate is M = 3.8 Mq Myr~^. 
At later times, the amount of transported dust varies around ^ 10% of the 
total dust mass. These variations are mainly due to dust particles born inside 
the corotation circle in the interarm region, where they are trapped into peri- 
odic orbits as shown by solid and dotted lines in Fig.5. We conclude that only 
a (small) fraction of dust particles (approximately 10% by mass) injected in- 
stantaneously inside the corotation circle can move to the outer regions, while 
the rest is confined in the inner galaxy. As a consequence, the radial transport 
of an instantaneously injected portion of dust terminates after 1 Gyr. 

The spiral stellar density wave considered so far has a pitch angle oi i — 25°. 
This is typical for Sc galaxies (Kennicutt, 1981). The efficiency of outward 
radial transport of dust may depend on the openness of stellar spiral arms. 
To test this conjecture, we explored a series of models with different pitch 
angles. The dashed line in Fig. 8 shows the the dust mass outside corotation 
as a percentage of the total dust mass for the stellar spiral with a pitch angle 
i — 15°. Roughly 8% of the injected dust mass is pushed outside corotation, 
as compared to ~ 10% for i = 25°. This indicates that the efficiency of out- 
ward radial transport of dust drops by roughly 2%. We expect that the radial 
transport of dust by spiral stellar density waves may becomes inefficient for 
very tightly wound spirals with a pitch angle i ^ 5°. Our simulations have 
also shown that the efficiency of outward transport of dust is weakly sensitive 
to the angular velocity Qgp of stellar spirals. For instance, a stellar spiral with 
Qsp = 21 km s^"*^ kpc^-*^ can transport roughly 1% less dust (sec the dotted 
line in Fig. 8) than a spiral with Q^p = 24 km s~^ kpc~^. We note that the 
corotation radius for figp = 21 km s~^ kpc~^ is at 8.8 kpc. 

The solid line in Fig. 8b shows the gas mass inside corotation as a percentage of 
the total initial gas mass. ^ It is evident that the gas (contrary to the dust) does 



^ Due to the imposed outflow outer boundary condition, we cannot reliably cal- 
culate the gas mass outside corotation because its value can be affected by the 
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Fig. 9. The ratio of the frictional timescale Tc to the dynamical timescale Tdyn as a 
function of radius at t = 1340 Myr. The solid and dashed lines show Tc/r^yn for the 
dynamical timescale defined as T^yn = ^^'^ "^dyn = l^d^ ~ ^sp^l' respectively. 
The detailed explanation is given in the text. 

not exhibit a noticeable radial migration. A very weak inward radial transport 
(~ 1% by mass) is seen for flgp = 24 km s~^ kpc~\ and it becomes negligible 
for Qgp — 21 km kpc~^. We attribute this difference between the gas and 
dust radial transport rates to the restoring force of the pressure gradients 
which are always present in the warm (Tg ~ 10^ K) gas. The warm gas resists 
radial redistribution by creating additional pressure gradients which act to 
restore the initial centrifugally balanced configuration. The dust, however, is 
cold and it can easily flow in the radial direction provided that the coupling 
between dust and gas is not very strong. 

To determine the main cause of the radial transport of dust, we compute 
the quantity Tc/rdyn (averaged over 0.5 kpc radial annuli). This quantity as 
a function of radius is shown by the solid line in Fig. 9 at the time when 
the gas spiral is fully developed. The coUisional timescale Tc of energy and 
momentum exchange between the dust and gas particles is defined in Eq. (4). 
The local dynamical timescale is defined as Tdyn = = r/vd,4„ where fd,,/, is 
the azimuthal component of dust velocity v^- If Tc <^ Tdyn the dust is expected 
to be strongly frictionally coupled to the gas and to essentially track the gas 

gas that leaves the computational area. The dust, however, never reaches the outer 
boundary. 
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motion. Figure 9 shows that Tc is a factor of 0.01 — 0.1 smaller than Tdyn, which 
implies a strong frictional coupling betwen dust and gas. However, judging 
from the bottom panels of Fig. 7, it would appear that the dust does not follow 
exactly the gas motion in the region near corotation. The reason is that the 
dust particles have natural resonant frequencies near corotation (remember 
that corotation determines the position of resonance in gas/dust orbits). If 
the gravitational field generated by spiral structure perturbs a dust particle's 
orbit at or near its resonant frequency, then the responce of the orbit will be 
large, even when the perturbing field is weak and/or the frictional couphng 
is strong. To investigate this effect, we re-define the local dynamical time as 
Tdyn = l^d^ ~^sp^l plot the corresponding quantity Tc/rdyn in Fig. 9 by the 
dashed line. It is now evident that the frictional coupling between dust and 
gas is much weaker near corotation. This is consistent with the trajectories 
of dust particles found in Sect. 4 - only the dust particles that are born near 
corotation show substantial radial migration. 



Another source of uncertainty in our simulations is the friction coefficient A, 
which we have adopted in the form of Eq. (4) to simplify the computations. 
As is pointed out in Sect. 2.1, the total stellar+gas surface density Eg + Eg 
(instead of simply gas surface density. Eg) should appear in the definition of 
the gas vertical scale height 



= 7rG(E;+Ej' ^^^^ 



which may considerably increase the value of A. The uncertainity in the value 
of Eg (we do not explicitly follow the evolution of the stellar disk) works against 
this re-definition of Zq in our simulations. In order to estimate the effect that 
the increased surface density of gravitating matter in Eq. (14) may have on the 
radial transport of dust particles, we notice that in the region of interest near 
corotation the ratio of stellar-to-gas surface densities is usually in the 1.5 — 15 
range. Therefore we performed a test run in which A was defined by Eq. (4), 
but was later multiplied by a factor of 10. The efficiency of outward radial 
transport of dust increased by less than 1%. In another test, A is decreased by 
a factor of 10, and the efficiency of outward radial transport of dust is dropped 
by approximately 1%. We conclude that an order of magnitude variation in A 
has a minor effect on the efficiency of radial transport of dust. This reinforces 
a previous conclusion: the resonant region near corotation, where ^ Qgp, 
and the collisional drag force is dynamically less important than the spiral 
gravitational field, determines the radial migration of dust. 



19 



6 Dust destruction 



In this section we study the effect of dust destruction by SNe on the efficiency 
of radial transport of dust in spiral galaxies. The dust grain destruction by 
SNe is the most important mechanism for cycling the dust back to the gas 
phase. The dust destruction rate can be determined as (Dwek, 1998) 

-Dd(r, 0, t) = mdest('^, 0, t) u{r, 0, t). (15) 



Here, ^{r, 0, t) is the supernova rate per unit area and mdest('", 0, t) is the mass 
of dust that is destroyed by a single SN expanding at location (r, 0) and time t. 
Only SN remnants with velocities exceeding 100 km s^^ can efficiently destroy 
dust grains. Therefore, the value of mdest(^) 4>j t) is usually estimated from the 
relation mdest/^d = eM^o^/Mg, where M^o" = 3.8 x 10=^ Mq (Lisenfeld and 
Ferrara, 1998) is the mass of a remnant accelerated to velocities > 100 km s~^ 
by the blast. This relation simply states that the fraction of destroyed dust 
(by mass) within a computational cell is directly proportional to the fraction 
of accelerated remnant mass. The efficiency of dust destruction is e ~ 0.1 
(McKee, 1989). 

The SN rate per unit area is defined as 

mup 

/ SsFR(r,0,t -r(m)) m-" dm 
r^{r,cl^,t) = '^ , (16) 

miow 



where Ssfr is the star formation rate per unit area, a is the slope of initial 
mass function (IMF), r(m) is the lifetime of a star of mass m, miow and mup 
are the lower/upper cutoff masses, respectively. In practice, we approximate 
the integral in Eq. (16) by the sum with a mass discritization dm — 1 M©. 

Observations of both normal and starburst disk galaxies suggest that on the 
scales of a few kiloparsecs star formation may be represented by a Schmidt 
law (Kennicutt, 1998) 



We modify the star formation law (17) by assuming that star formation is 
suppressed when Eg < min{i:^9°'^'^ ,11^}^^'). 

Often, the instantaneous recycling approximation is used to solve Eq. (15), 
which neglects the time delay that enters Eq. (15) through the SN rate v. 
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However, wc follow a more accurate approach and form stellar clusters each 
Myr according to the star formation law (17). Stellar clusters are assigned po- 
sitions and velocities drawn from the parent gas and are evolved as coUisionless 
particles in the combined gravitational potential of the halo and stellar disk. 
The local self-gravity among stellar clusters is neglected. Since each stellar 
cluster carries information on the star formation rate at the time of its birth, 
the supernova rate z/(r, and the dust destruction rate Dd{r,(f),t) in each 
computational cell (r, 0) at a time t can be easily obtained. This approach 
was successfuUy applied by Vorobyov (2003) to model the Ha luminosity in 
ring galaxies. We use the Salpeter IMF with a — 2.35 and lower/upper cutoff 
masses miow = 0.5 M© and mup = 40 Mq, respectively. 

We modify the hydrodynamic equations (l)-(3) of our gas+dust system to take 
into account star formation and dust destruction. The continuity equations for 
the gas and dust become 



9Eg 
dt 

dt 



+ V-{^,v,)=S,-D„ (18) 
+ V-(Edt;d) = -^d. (19) 



The resulting momentum equations are 



^ SgVg + (vg • V)EgVg = -EgV$si,s2,h - VPg + Ed/ + 

v,{S,-D,), (20) 

^ Eava + (Vd • V)Edt;d = -Ed V$si,s2,h - Ed/ - v^D^. (21) 



We have rewritten equations of motion in the form of momentum equations for 
E V which are appropriate for the system with sources and sinks. Star forma- 
tion depletes the gas reservoir of our model galaxy and is taken into account 
by the sink term Dg — EgpR in Eqs. (18) and (20). The gas is returned to the 
system with supernova explosions and quiet mass loss of intermediate and low 
mass stars. Since we do not use a multiphase description of the interstellar 
medium, the ejected hot gas is directly transformed into a warm (T ~ 10"^ K) 
phase. The rate of gas ejection per unit area is determined by 



^ miow 



/ g{m) EsFR(r, 0, t - r(m)) m dm 



mup 
miow 



(22) 
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where g{m) denotes the gas mass ejected by a star of mass m (Koppen and Ari- 
moto, 1991). We found it computationally prohibitive to compute accurately 
(i.e. by taking into account the time delay r(m) in Eq. [22]) the gas ejection 
rate by stellar clusters on a time scale of interest (a few Gyr). Therefore, for 
the stellar clusters older than 100 Myr, the instantaneous recycling approxima- 
tion is assumed by neglecting the time delay in Eq. (22). These stellar chisters 
instantaneously release all gas that can be produced by the 0.5 — 5 Mq stars. 
Further dynamical evolution of such clusters is not computed. 

We study the effect of dust destruction on the radial transport of dust using 
a prototype model in which the spiral density wave has a pitch angle i = 25° 
and angular velocity fl = 24 km s~^ kpc~^. This model was studied in detail 
in Sect. 5 where we neglected the effect of dust destruction. It was shown there 
that about 10% of dust by mass injected instantaneously inside the corota- 
tion circle can be transported outside this region by the gravitational drag 
associated with spiral density waves. We reproduce this result in Fig. 10a by 
the dashed line, and the solid line gives the dust mass outside corotation as a 
percentage of the initial total dust mass when dust destruction is taken into 
account. As is seen, supernova explosions have a minor effect on the radial 
transport of dust, decreasing the relative amount of dust outside corotation 
by only 2% when compared to the model without dust destruction. It is worth 
noting that inside corotation the destruction of dust by supernova explosions 
is efficient. Indeed, the dashed-dotted and dotted lines in Fig. 10b show the 
dust mass inside corotation as a percentage of the initial total dust mass when 
dust destruction is taken and not taken into account, respectively. Approxi- 
mately 15% of the instantaneously injected dust mass is left inside corotation 
after 2 Gyr. Given that another 10% has been transported outside corotation, 
SNe have destroyed about 75% of the injected dust mass during 2 Gyr. In 
Sect. 5 we argued that the ultimate fate of dust grains would depend on their 
place of birth. The simulations discussed in this section reinforce our earlier 
conclusion. In our model dust grains are formed inside the corotation circle at 
~ 8 kpc. The dust grains that are born in the inner 6.0 kpc circle (by im- 
plication they constitute the majority of the total dust mass) move on elliptic 
orbits periodically crossing the spiral arms and will be inevitably destroyed 
by supernova explosions after several revolution periods. However, the dust 
grains that are born in the 6 kpc <^ r < 8 kpc annulus and sufficiently away 
from spiral arms (by implication they constitute the minority of the total dust 
mass) manage to survive the destructive effect of supernova explosions during 
many revolution periods of our model galaxy. These dust grains may move on 
orbits similar to those shown in Fig. 5 by the solid and dotted lines, circulating 
between the spiral arms (but never crossing them) and periodically traveling 
in and out of the corotation circle. We conclude that a (small) portion of dust 
grains may survive the destructive effect of SNe during many (> 10) revolution 
periods of a spiral galaxy. 
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Fig. 10. a) Dust mass outside and b) inside corotation as a percentage of the total 
initial dust mass. The solid and dashed lines show the relative dust mass outside 
corotation in the models with and without dust destruction, respectively. The dotted 
and dotted-dashed lines give the relative dust mass inside corotation in the models 
with and without dust destruction, respectively. 
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7 Summciry 



In this paper we study the outward radial transport of interstellar dust grains 
which is connected with the perturbation of circular motions caused by spiral 
stellar density waves in disk galaxies. We consider a model spiral galaxy in 
which star formation (and, by implication, dust formation) is localized within 
the corotation circle. We find that the combined action of the gravity force 
of spiral stellar density waves and the drag force between the gas and dust 
components imparts to dust particles a radial velocity component, and thus 
can evacuate dust outward. We conclude that: 

1. Dust grains that are formed inside corotation can be transported to a dis- 
tance that exceeds the corotation radius by roughly 25%. In particular, if coro- 
tation is located at Vcr = 8.0 kpc, dust grains formed inside the 6 kpc <^ r < 
8.0 kpc annulus can travel to distances as large as 10 kpc in roughly 1.0 Gyr. 
This is approximately an order of magnitude faster that can be provided by 
the interstellar turbulence. 

2. A fraction of dust grains can be trapped into the interarm region, and thus 
can escape a hostile environment from supernova explosions acting mostly in 
the spiral arms. These grains form diffuse spiral arms which extend 4-5 kpc 
from the sites of active star formation. Dust particles that cross the spiral 
arms during their drift outward can survive if they do not meet strong shocks 
from SNe while passing through the arm. 

3. In the model that considers instantaneous dust injection, about 10% of dust 
by mass is evacuated outside the corotation circle in roughly 1.0 Gyr. One can 

assume that continuous injection will result in the same percentage of evacu- 
ated dust. The efficiency of outward radial dust transport is weakly sensitive 
to the rotation speed of stellar spirals, but it may be considerably reduced for 
very tightly wound spirals with a pitch angle i ^ 5°. Dust destruction by su- 
pernova explosions has only a minor effect on the efficiency of outward radial 
transport of dust by spiral density waves, which decreases the relative amount 
of transported dust mass by ~ 2%. A (small) portion of dust grains may sur- 
vive the destructive effect of SNe during many (> 10) revolution periods of a 
spiral galaxy. 

4. Our modeling shows that the dust-to-gas ratio in the region near corotation 
(±2 kpc) can vary substantially over the scales comparable to the thickness 
of spiral arms, ~ 1 kpc. Possible observations of such variations would be 
suggestive of the proposed mechanism of radial dust transport. 

In our model we neglected the fact that charged dust grains are strongly 
coupled to the magnetic field - the gyration frequency for a typical grain 
with radius a ~ O.l/zm and grain charge Z ~ 100 is of cob ~ 10~^° Hz. We 
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expect that in more realistic models, where the magnetic field is permanently 
regenerated by a turbulent dynamo, dust particles will be involved in both 
the turbulent difi:usive motions and regular flow, and the overall picture will 
remain qualitatively similar to that described above. We leave this issue for 
further study. 

We should mention that the interaction of gas and solid particles through the 
drag force and the associated radial flow of solid particles have been already 
extensively studied in the context of migration of planetisimals in protoplan- 
etary disks. For example, Weidenschilling (1977) has noticed that in proto- 
planetary disks the planetesimals should migrate in the direction of increasing 
pressure (i.e. inward, for an axisymmetric disk with a declining radial density 
proflle). This is due to the fact that the gas in protoplanetary disks generally 
orbits the central protostar at the sub-Keplerian velocity due to a substantial 
contribution of pressure gradients in the gas rotation curve. The cold planetes- 
imals (which by implication have no pressure) are then forced to flow inward 
if the frictional orbital coupling between the gas and planetesimals is not too 
strong. This mechanism however may not work on galactic scales because 
the contribution of the pressure gradients to the gas rotation velocity may 
be negligible (see Sect. 2.2). In fact, the pressure gradients may even assist 
the gravity in some disk galaxies with a ring-like radial distribution of gas. 
Another interesting aspect of gas-planetesimal interaction in self-gravitating, 
protoplanetary disks was investigated by Rice et al. (2004). They considered 
disks with a flocculent spiral structure and demonstrated that the drag force 
rather than gravitational field of spiral arms plays a dominant role in the dy- 
namics of planetesimals of intermediate size, forcing them to concentrate in 
the gas spiral arms. The non- axisymmetric gravitational field of grand-design 
spirals (considered in the present paper) is by implication much stronger than 
that of flocculent spirals. As a consequence, the trajectories of dust particles 
can be considerably influenced by the spiral gravitatioanl field (Fig. 5). We 
believe that the resonant action of spiral gravitational field on dust particles 
(rather than the drag force between the gas and dust) is the main driving 
force for the radial outward transport of dust in grand-design spiral galaxies. 
This conjecture is justified by the fact that the radial fiow of dust takes place 
near the corotation resonance where a weaker frictional coupling between the 
dust and gas occurs (see Fig. 9). 
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Abstract 

Motivated by recent observations which detect dust at large galactocentric distances 
in the disks of spiral galaxies, we propose a mechanism of outward radial transport 
of dust by spiral stellar density waves. We consider spiral galaxies in which most of 
dust formation is localized inside the corotation radius. We show that in the disks 
of such spiral galaxies, the dust grains can travel over radial distances that exceed 
the corotation radius by roughly 25%. A fraction of the dust grains can be trapped 
on kidney-shaped stable orbits between the stellar spiral arms and thus can escape 
the destructive effect of supernova explosions. These grains form diffuse dusty spiral 
arms, which stretch 4-5 kpc from the sites of active star formation. About 10% of 
dust by mass injected inside corotation, can be transported over radial distances 
3 — 4 kpc during 1.0 Gyr. This is roughly an order of magnitude more efficient 
than can be provided by the turbulent motions. 

Key words: galaxies: spiral, ISM: dust 
PACS: 98.58.Ca, 98.52.NR 



1 Introduction 

Interstellar dust has been recently recognized to spread over large distances 
from their source, and thus can be considered as a tracer of dynamical pro- 
cesses responsible for the circulation of material in disk galaxies. Observa- 
tions of edge-on galaxies reveal a great amount of dust far outside the galactic 
planes: in NGC 891 and NGC 4212 dusty clouds extend up to z ~ 2 kpc above 
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the plane, and contain in total as much as ~ 10^ Mq (Howk and Savage, 1997, 
1999; Rossa et al., 2004, Thompson et al., 2004). Several mechanisms can con- 
tribute to transport of dust in the vertical direction: convective motions or 
bore flows (Gomez and Cox, 2002) associated with spiral density waves, chim- 
neys produced by multiple supernova explosions (Norman and Ikeuchi, 1989), 
and radiation pressure (Ferrara et al., 1991; Shustov and Vibe, 1995; Dettmar 
et al., 2005). 



The situation is less clear when the radial distribution of dust is considered. 
Observations of spiral galaxies in the far infrared band show that dust can ex- 
tend in the radial direction far beyond the stellar disks. Neininger et al. (1996) 
found that the radial distribution of 1.2 mm dust emission from NGC 4565 is 
two times wider than CO emission. ISO observations of several spiral galaxies 
at 200 nm (Trewhella et al., 2000; Bianchi et al., 2000) revealed the presence 
of relatively cold dust outside the stellar disks with a scale length i?^ ;^ 1.5 i?* 
{Rd and are the radial scales of dust and stars, respectively). A similar 
conclusion was made by Xilouris et al. (1999) for four spirals observed in the 
optical and near infrared bands. Bianchi et al. (2000) found that the best flt 
of the 200 nm emission for NGC 6946 gives Rd ^ 3 R^,. MHD turbulence can, 
in principle, spread dust outward in radial direction (Cho et al., 2003), how- 
ever, the corresponding time for galactic scales seems to be too long. Indeed, 
the coefficient of turbulent diffusion D ~ {Lv)/3 is of the order of 10^^ cm^ 
s"^ for typical interstellar values L ~ 100 pc and ^; ~ 10 km s^. Therefore, 



dust particles can diffuse by turbulent motions over scales ^(Ai22) > 3 kpc 
in a characteristic time t > 16 Gyr. In these conditions, only regular hydrody- 
namic motions can explain the presence of dust outside stellar disks in spiral 
galaxies. 



In this paper we explore the possibility that the radial transport of dust in 
spiral galaxies is driven by the hydrodynamic motions associated with spiral 
stellar density waves. In Sect. 2 a set of model hydrodynamic equations for 
dusty interstellar medium is described and the equilibrium parameters are 
flxed. In Sect. 3 the generation of spiral pattern in the stellar disk is formulated. 
Sect. 4 describes trajectories of single dust particles and Sect. 5 describes the 
dynamics of dust fluid as a whole including formation of dusty spiral arms. 
Sect. 6 summarizes our results. 
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2 Model description 



2.1 Basic equations 



Our model disk is a two-fluid system composed of dust and gas, which evolves 
in the external potential of the stellar disk and dark matter halo. We assume 
that the dust and gas are coupled by friction. The gas disk is isothermal at T = 
10^ K and the dust pressure is negligible. We use the thin disk approximation, 
i.e. all quantities are vertically integrated. The system is described by the basic 
hydrodynamic equations, where we denote the dust and gas by the subscripts 
d and g, respectively. 



The continuity equation for the gas and dust is 

+ V-(Sg,dt;g,d)=0, (1) 



dt 

where Eg^d are the gas and dust surface densities and 'yd,g are their velocities. 
The equations of motion are 



^ + K • V)t,g = - V$.,,,,H - ^ + (2) 

^ + (vd • V)i;d = - V$si,s2,h - /, (3) 

where, $si,$s2j and $h are the contributions to the total gravitational po- 
tential by the axisymmetric stellar disk, non-axisymmetric spiral density per- 
turbation, and dark matter halo, respectively. The friction force is defined as 
/ = A{vd - Vg) (see Draine and Salpeter, 1977; Noh et al., 1991). 

We take the friction coefficient A to be a simple function of the disk conditions. 
We assume that the thickness of the dust disk is comparable to that of the 
gas disk. Then the coUisional timescale Tc of energy and momentum exchange 
between the gas and dust particles may be used as an estimate of A~^, 

. 1 (TcngVthmg 

A^T^ = , (4) 

TOd 

where Ug is the number density of gas particles, TOg and TOd are the masses of 



gas and dust particles, and Vth = y3RT/ fi is the gas velocity dispersion; Tc is 
estimated for the subsonic relative motion of neutral dust grains and gas. Since 
we use the thin disk approximation, the gas volume density in Eq. (4) is ap- 
proximated assuming a local vertical hydrostatic equilibrium nigHg = T,g/zo, 
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where the gas disk thickness is defined as^ zq — 2Cg/(7rGS„), where Cg is the 
sound speed. The geometrical cross section is used as an estimate of the coUi- 
sional cross section ctc = vra^, where a is the average radius of a dust particle. 
In the present simulations, we assume ma = 10~^^ g and a — 10~^ cm. It 
is well known that interstellar dust grains are charged, and Coulomb inter- 
actions with charged ambient ions and electrons increase the friction force 
(Draine and Salpeter, 1979; Weingartner and Draine, 2001). However, as we 
are mostly concerned with the diffuse HI phase of the ISM (warm neutral 
medium, WNM, in standard nomenclature), this increase does not crucially 
change the resulting friction coefficient. Indeed, the Coulomb friction can be 
written as Fc = x(^)-^hi! where F^j is the friction in an HI environment, 
x{Z) — {az/ayiriAJ2xi{mi/'mp), az = Ze'^/kT, In A is the Coulomb loga- 

i 

rithm, Xi is the fractional ionization of i-th ions or electrons, is the corre- 
sponding mass, averaging is over charge distribution of dust grains (Weingart- 
ner and Draine, 2001). For typical conditions in the WNM; a; ~ 0.1, n ~ 0.3 
cm~^, and dust charge (Z) ~ 60 (Yan et al., 2004) we get x ~ 1. Therefore, 
we will restrict our consideration with the friction coefficient A as determined 
for the neutral dust. 



2.2 Equilibrium gas disk 



Initially, our model galaxy consists of a rotating gas disk balanced by the grav- 
ity of an axially symmetric stellar disk and spherically symmetric dark matter 
(DM) halo. The dust disk is introduced later when the non-axisymmetric part 
of the stellar gravitational potential is set. The rigid stellar disk is assumed to 
have a power-law radial density profile of the form 



27rG 



(5) 



where B'^ = 27rG'rg Ego. The gravity force of such a density distribution is given 
by Toomre (1963) 



dr 



^B^ 



+ r') 



-3/2 



(6) 



In the following, we use the central stellar density Ego = 1.2 x 10^ Mq pc ^ 
and rs — 3 kpc, which gives us a total stellar mass of Mgt — 7 x 10^^ Mq. 

^ Since the dominant component of gravity in the vertical direction is the stellar 
disk, the total surface density (Stot = + ^s) should appear in the expression for 
zq. This may increase the friction coefficient A, especially in the inner regions. We 
discuss this issue in Sect. 5 
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The density distribution of the rigid DM halo is assumed to be that of a 
modified isothermal sphere (Binney and Tremaine, 1987) 

(1 + r/rhY 



where the central density pho and the characteristic scale length are given 
by Mac Low and Ferrara (1999) and Silich and Tenorio-Tagle (2001): 

Pho = 6.3 X 10i« (^] h~"^ Mo kpc-=^ (8) 



, /M- \l/2 

rh = 0.89 X 10-' ( ^ ) h^l^ kpc. (9) 



Here, h is the Hubble constant in units of 100 km s^^ Mpc~^ and Mh — 
10^^ Mq is the total DM halo mass. We note that the DM halo mass within 
the computational domain r < 20 kpc is Mh(r < 20 kpc) = 5 x 10^" Mq. We 
adopt h = 0.65 throughout the paper. The gravitational force of spherically 
symmetric DM halo can be expressed as 



dr 



AnGphoTh [r/rh - arctan(r/rh)] (^^^ . (10) 



The gas disk has an exponentially declining density profile with the central 
surface density Ego = 30 M© pc~^, and radial scale length rg = 9 kpc. The 
total mass of the gas disk within the computational domain (r = 20 kpc) is 
Mg = 1.0 X 10^" Mq. The gas disk contains only a small fraction (~ 10%) of 
the total mass in the computational domain. Hence to a first approximation 
we can assume that the gas moves in the potential field of the stars and DM 
halo. 

Once the density profile of the gas disk is fixed, its initial rotation curve 
(RC) is obtained by solving Eq. (2), where the friction force / is set to zero 
and the time variations are neglected. The resulting initial RC is plotted in 
Fig. 1, where it reaches a maximum circular velocity of 200 km s~^ at r fa 
3 kpc. Beyond this maximum, the rotational velocity undergoes a gradual 
decline. Because the gas surface density has a slowly declining radial profile 
(rg = 9 kpc) and the gas is assumed to be isothermal, the contribution of the 
pressure gradient in the gas rotation curve is less than 1% everywhere in the 
disk. 



5 



250 




I ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' ' ' ' ' ' 1 

5 10 15 20 

Radial distance (kpc) 



Fig. 1. The initial rotation curve of the gas disk. 
2.3 Computational techniques 

A set of hydrodynamic equations (l)-(3) in polar coordinates (r, 4>) is solved 
using the method of finite-differencing with a time-explicit, operator-split so- 
lution procedure (ZEUS-2D) described in detail in Stone and Norman (1992). 
We use a resolution of 200 x 200 grid points in a polar computational domain 
with a radius of 20 kpc. We set a free boundary condition at the outer bound- 
ary, i.e. the gas and dust are allowed to flow out freely from the computational 
domain. The outflowing gas/dust is assumed to be lost by the parent galaxy. 



3 Generation of spiral pattern 

Once the axisymmetric equilibrium gas disk is constructed, we introduce the 
non-axisymmetric spiral gravitational potential of the stellar disk $s2- The 
latter is deflned in the form of a running spiral density wave (Lin et al., 1969) 

*s2 = -C{r) cos [m(cot(i) ln(r/rsp) + (p - Qspt)] ■ (11) 
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Fig. 2. Behaviour of n±«;/2, where k is the epicychc frequency and Q is the angular 
velocity of gas. The angular velocity of spiral pattern figp is shown by the solid line. 

Here, C(r) is the amplitude of the spiral gravitational potential, m is the 
number of spiral arms, i is the pitch angle, rgp is the characteristic winding 
length of spiral arms, and Ogp is the angular velocity of spiral pattern. In the 
following we adopt m = 2, i = 25°, Vgp = 6 kpc, flgp = 24 km kpc~^. 
This choice of Ogp places the corotation of the gas disk at ~ 8.0 kpc and 
the outer Lindblad resonance at fa 12.0 kpc. The inner Lindblad resonance 
is absent. The behavior of gas angular velocity Q and Q ± k,/2 is shown in 
Fig. 2. According to Kennicutt (1981), a pitch angle oi i = 25° is typical for 
Sc galaxies. The efficiency of radial transport of dust for different pitch angles 
is discussed in Sect. 5. 

The logarithmic spiral defined by Eq. (11) is quickly winding up as r — > 0. 

This may create large gradients of gravitational potential near the origin and 
lead to serious numerical difficulties. For instance, if the amplitude C(r) is 
independent of r, than the ratio of maximum non-axisymmetric perturbing 
force to the total axisymmetric gravity force /3{r) — max |V$s2|/(5$si/5r + 
d^h/dr) sharply increases as r — )■ 0. This is illustrated in Fig. 3, where /3{r) 
is plotted by the dotted line for C = 0.005 (in dimensionless units). The 
apparent growth of (3 with radius for r > 3 kpc is also unrealistic since the 
spiral arms are expected to vanish at larger radii. In view of these difficulties, 
we use a radially dependent amplitude of the form C(r) = Co (r )"(''). Here, 
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Fig. 3. The ratio of the maximum, non-axisymmetric perturbing force 
|V$s2| = [{d^s2/dr)'^ + (r~-^9$s2/9<?!))^]^/^ to the total a:xisymmetric gravity force 
{d^si/ df + d^h/dr) as a function of radial distance. The solid and dashed lines are 
explained in greater detail in the text. 

Co(r) is a linear function of r which has a value of at r = kpc (ensuring 
that |V$s2| diminishes as r — )■ 0) and attains its maximum value of 0.0045 
(in dimensionless units) at r = 20 kpc. The exponent a{r) decreases linearly 
with radius from a = 2 at r=0 kpc to a = —2.1 at r=20 kpc. The resulting 
profile of /3(r) is shown in Fig. 3 by the solid line. Since the non-axisymmetric 
gravity force scales as |V$s2| oc l/r-^~"*^''\ the ratio P increases with radius 
at r < 4 kpc, and decreases at r ;^ 4 kpc. At the position of corotation 
Tcr ~ 8 kpc, (3 drops by roughly a factor of 2 as compared to its maximum 
value at r ~ 4 kpc and continues to decrease at larger radii. We note that the 
maximum, non-axisymmetric perturbing force never exceeds 12% of the total 
axisymmetric gravity force. 

The non-axisymmetric part of the stellar gravitational potential $s2 is turned 
on slowly. Specifically, $s2 is multiplied by a function e{t), which has a value 

of at t=0 and linearly grows to its maximum value of 1.0 at t > 200 Myr. 
It takes a few hundred Myr for the gas disk to adjust to the spiral distortion 
and develop a spiral structure as shown in Fig. 4. The gas spiral pattern (as 
well as the stellar one) rotates counterclockwise. The position of corotation 
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Fig. 4. Gas density distribution att = 0.8 Gyr from the beginning of numerical sim- 
ulations. The corotation radius is shown by the circle. The scale bar is in Mq pc~^. 
The contour line sketches the region of supercritical gas density (see Sect. 5 for 
details). 

is shown by the circle. The strongest gas response to the underlying stellar 
spiral density wave is seen inside the corotation circle. However, a weak spiral 
structure can still be traced outside corotation at radii r = 10 — 15 kpc. 



4 Trajectories of dust pzirticles 

In this section we study the trajectories of single dust particles in the com- 
bined gravitational potential of stellar disk and dark matter halo including 
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Fig. 5. Polar plot (r, 6) showing trajectories of three test dust particles. The radial 
distance is in kiloparsecs and the azimuthal angle is in radians. Stellar spiral arms 
are located along azimuthal angles of and tt. The direction of motion of test 
particles in the frame of reference of stellar spiral arms is indicated by arrows. 

the effects of dust-gas friction. The equation of motion of the dust fluid (3) 
can be obviously reduced to a set of equations of motion for a single dust par- 
ticle and subsequently solved using the Runge-Kutta scheme. We introduce 
test dust particles at different radial distances from the galactic centre and 
at different azimuthal angles 6d with respect to the position of stellar spiral 
arms. In our simulations, the stellar spiral arms are modeled by the spiral 
gravitational potential rotating counterclockwise in the laboratory frame of 
reference. Therefore, the position of stellar spiral arms at any given time is 
determined as the position of minimum in the spiral gravitational potential 
defined by Eq. (11). All dust particles are injected att = 0.8 Gyr when the gas 
spiral structure is fully developed (see Fig. 4), and they are assigned velocities 
drawn from the parent gas. We followed the trajectories of test dust particles 
for 2 Gyr. 

The phase (ra, ^d) trajectories of three representative dust particles in the 
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Fig. 6. Schematic description of gravity forces (shown by the dashed arrows) acting 
on test dust particles from the stellar spiral arms (plotted by thick solid and dashed 
lines) . Velocities of dust particles in the frame of reference of stellar spiral arms are 
indicated by solid arrows. Two test particles inside corotation and four test particles 
outside corotation are shown. Arrows are not to scale. 

frame of reference of stellar spiral arms are shown in Fig. 5. The stellar spiral 
arms in Fig. 5 are always located along azimuthal angles — and tt. 
The dashed line shows the trajectory of a dust particle placed initially at 
rd=5.0 kpc in the arm (^d = 0; hereafter, particle 1), whereas the solid line 
gives the trajectory of a dust particle placed initially at rd=7.0 kpc in the 
interarm region (^d = 7r/4; hereafter, particle 2). Dust particles similar to 
particle 1 may be injected by supernova explosions which tend to occur near 
spiral arms. Particle 2 may be created by red giant stars. 

It is apparent that the trajectories of particles 1 and 2 are quite different. 
Particle 1 moves on an elliptical orbit in one direction (counterclockwise), re- 
maining always inside the corotation circle (vcr ~ 8.0 kpc). The trajectory 
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of particle 1 is fairly stable. It comes closer to the origin when traveling be- 
tween the arms (^d ~ 7r/2) and it passes through the spiral arms 4 times 
during 0.65 Gyr. Since dust particles are thought to be destroyed by super- 
nova explosions, such a frequent encounter with the spiral arms implies that 
particle 1 may not live longer than 0.5 Gyr (although it depends on the cov- 
ering factor of sufficiently strong shock waves, vstss ^ 100 km s~^, generated 
by supernova explosions (SNe) and propagating perpendicular to the spiral 
arms). Particle 2 has a stable, kidney-shaped trajectory, which shows an im- 
pressive radial migration of ~ 4 kpc, as opposed to only 1 kpc for particle 1. 
Interestingly, particle 2 is caught between the spiral arms, periodically passing 
through the corotation circle from the inner part of the disk to the outer part 
and vice versa. Since particle 2 never encounters the spiral arms during 2 Gyr, 
it may live significantly longer than particle 1. Dust particles that are born 
at ~ 7''/2 (i.e. exactly between the spiral arms) but are sufficiently away 
from the galactic center, move on much more compact orbits as shown by the 
dotted line in Fig. 5. 

To summarize, we find that all particles born at r < 6 kpc, i.e. ^ 2.0 kpc 
away from the position of corotation radius Tcv ~ 8.0 kpc, always stay inside 
corotation. This is irrespective of their azimuthal angles of birth. Trajectories 
of those dust particles resemble that of particle 1 shown in Fig. 5 by the dashed 
line. Such particles often encounter spiral arms and thus may be destroyed in 
one rotation period by shocks from SNe. The fraction of surviving dust grains is 
determined by the probability for them to meet a sufficiently strong {vg ^ 100 
km s~^) destroying shock wave inside the arm. On the other hand, particles 
injected near corotation at 6 kpc <^ r < 8.0 kpc (and sufficiently away from 
spiral arms) spend most of their lifetime in the interarm region. Therefore, 
most of these dust particles can survive over the entire migration period in 
the outer galactic regions. 

In the absence of spiral arms, dust particles move on circular orbits deter- 
mined by the axisymmetric global gravitational potential of a galaxy and the 
drag force from the gas disk. When the non-axisymmetric spiral gravitational 
potential is introduced, dust particles move on various non-circular orbits as 
shown in Fig. 5. In order to understand the origin of such non-circular orbits, 
we consider the influence of a non-axisymmetric spiral gravity force on the 
motion of dust particles. The gravity force of spiral arms acting on moving 
particles is schematically plotted in Fig. 6 by the dashed arrows. Six parti- 
cles are presented; two particles inside the corotation circle and four particles 
outside the corotation circle. The direction of motion of dust particles in the 
frame of reference of the spiral pattern is shown by the solid arrows. Two 
spiral arms are schematically shown by the solid and dashed lines. A moving 
dust particle may see the nearest spiral arm as either convex or concave. If 
the dust particle is closer to the concave spiral, the latter acts to increase 
(on average) the radial distance of the particle due to its gravitational drag. 
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Fig. 7. Gray-scale images of the gas distribution (the upper panels), dust distribution 
(the middle panels) and the dust-to-gas ratio (lower panels). The elapsed time is 
indicated in each panel. The position of corotation of the gas disk is outlined by 
the circle. The scale bars in the upper and middle panels give the column density 
in Mq pc~^ 
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On the contrary, if the dust particle is closer to the convex spiral, the latter 
tends to decrease the radial distance of the particle. If we now recall that the 
trailing spiral arms in Fig. 5 are located along the azimuthal angles and 
TT and rotate counterclockwise, then the effect of gravitational drag from the 
spiral arms becomes clearly evident in the trajectories of dust particles. For 
instance, when dust particles are in the first quadrant, they are attracted by 
the convex spiral located at = 0. As a consequence, the radial position of 
these particles decreases. When dust particles enter the second quadrant, their 
motion is controlled by the concave spiral located at ^sp = and their radial 
position increases. Therefore, the trajectories of dust grains may be sensitive 
to the pitch angle of stellar spiral arms. The maximum radial migrations of 
dust grains are expected for a pitch angle of 45°. 

The action of the friction force is more difficult to assess. The friction between 
dust and gas may help reverse the direction of motion of a dust particle (with 
respect to the spiral arms) when it passes through corotation. Indeed, when 
the dust particle is inside the corotation circle, its velocity is larger than that 
of spiral arms (the same is true for the gas in general) . After passing through 
corotation, the dust particle enters the region which is characterized by a 
smaller velocity than that of spiral arms. Therefore, we may expect that the 
friction between gas and dust would act to equalize the velocities of dust 
and gas, and thus decelerate the dust particles. This would eventually reverse 
their direction of motion with respect to the spiral arms as seen in Fig. 5 for 
particle 2. 



5 Response of dust to the stellzir spiral density wave 



In this section we study how dust responds to the imposed spiral disturbance. 
The dominant sources of dust formation in disk galaxies are stellar winds from 
late-type giant and supergiant stars, supernova explosions, and growth of pre- 
existing grains in dense n ^ 10^ cm~^ molecular clouds (Dwek, 1998). All 
of these sources are mainly located inside the corotation circle. Indeed, the 
strongest gas response to the stellar spiral density wave in Fig. 4 is obviously 
seen inside corotation. We compute the critical density for star formation as 
defined by the Toomre criterion (Kennicutt, 1989) 



vToomre _ ^-^ Cg K, 

^crit — _ ri ' ^ ^ 
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where k is the epicychc frequency, and the shear criterion (Martin and Ken- 
nicutt, 2001) 

V shear _ 2.5Ash f s /-, o\ 

^^^^ 



where = —0.5 r dO/dr is the local shear rate (the Oort constant). The 
contour line in Fig. 4 traces the region of supercritical density with Eg > 
min(E(JJ;°j°™'''^, Ej.^j.^^'^'') which is evidently localized inside corotation. Although 
the gas density outside corotation may occasionally become supercritical (not 
shown in Fig. 4) especially along the outer diffuse gas spirals, a lower gas 
density and higher temperature as compared to the inner region can work 
against star formation in that region. We thus conclude that most of star 
formation in our model galaxy is localized within the corotation circle. The 
radial scale length of our old stellar disk is rg = 3 kpc which also implies that 
most of low to intermediate mass stars that are capable of producing dust are 
confined within the corotation circle. 

We instantaneously inject 3.4 x 10^ Mq of dust mass inside the corotation 
circle (r < 7.8 kpc) and then determine the fraction of dust that can move 
outside corotation. In the present simulations, the dust particles are injected 
at t = Myr when a spiral structure in the gas disk has not yet developed. 
We also tried to introduce dust at t=800 Myr (when the gas spirals are fully 
developed) and confirmed that the time of dust injection has little infiuence 
on dust dynamics. The spatial distribution of injected dust is such that the 
dust-to-gas ratio 7 = 0.01 inside corotation, and 7 = outside corotation. The 
velocity of injected dust is equal to that of the gas. A conservative estimate 
of dust-to-gas mass ratio in the solar neighborhood gives an average value of 
7 = 0.006 ranging from 0.002 up to 0.04 in H2 regions (Spitzer, 1978). 

The upper and middle panels in Fig. 7 show two snapshots of the evolution of 
the gas and dust surface densities, respectively, whereas the lower panels give 
the dust-to-gas mass ratio. The numbers in each panel indicate the elapsed 
time since the beginning of the simulation. The circle indicates the position 
of corotation of the gas disk. A strong spiral response to the stellar density 
wave develops in the dust disk after approximately 3 revolution periods of our 
model galaxy. The dust density closely follows that of gas in the inner regions 
of our model galaxy. Small radial and azimuthal variations of the dust-to-gas 
ratio around the initial value 0.01 are indicative of a dominant role of frictional 
force in the dynamics of dust in the inner 5 kpc. This is in agreement with the 
results of Sect. 4 where we have found that dust particles born in the inner 
5 kpc move on fairly stable elliptical orbits around the galactic center. At 
larger radii 6 kpc ^ r < 8.0 kpc, however, significant variations of 7 become 
evident. There are regions near corotation that are completely devoid of dust, 
even though the initial dust-to-gas ratio is 0.01. A noticeable portion of dust is 
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Fig. 8. a) Dust mass ouside corotation as a percentage of the total dust mass. 
The sohd and dashed lines correspond to the models with pitch angles of stel- 
lar spiral density wave i = 25° and i = 15°, respectively. The dotted line gives 
the dust mass outside corotation for a lower rotation frequency of stellar spirals 
rigp = 21 km s~^ kpc~^. b) Gas mass inside corotation as a percentage of the initial 
total gas mass. The solid, dashed, and dotted lines correspond to the same models 
as in the upper panel. The difference in the values of the gas fraction inside corota- 
tion for ^Igp = 21 km s^^ kpc^^ and = 24 km s~^ kpc~^ is due to the different 
positions of the corotation radius in the two cases. 
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observed outside corotation. This dust tends to form into diffuse outer spiral 
arms, which appear to be a smooth extension of the inner dust spiral arms and 
can be traced up to r = 10 kpc. The dust-to-gas ratio in the outer dust arms is 
not dissimilar to that of the inner r < 5 kpc region. We thus conclude that the 
perturbing force of stellar spirals starts to play an important role in the dust 
dynamics near corotation resulting in an outward radial transport of dust. If 
we now identify the regions of supercritical gas density by contour lines in the 
upper panels of Fig. 7, than it becomes evident that the dust component can 
be found as far as 4-5 kpc from the sites of active star formation. 

The solid line in Fig. 8a gives the dust mass outside corotation (as a percentage 
of the total dust mass) versus time elapsed since the injection of dust. Approx- 
imately 10% of the initial injected dust mass is transported outside corotation 
during 1 Gyr. This implies that the dust transport rate is M = 3.8 Mq Myr^^. 
At later times, the amount of transported dust varies around ^ 10% of the 
total dust mass. These variations are mainly due to dust particles born inside 
the corotation circle in the interarm region, where they are trapped into peri- 
odic orbits as shown by solid and dotted lines in Fig. 5. We conclude that only 
a (small) fraction of dust particles (approximately 10% by mass) injected in- 
stantaneously inside the corotation circle can move to the outer regions, while 
the rest is confined in the inner galaxy. As a consequence, the radial transport 
of an instantaneously injected portion of dust terminates after ^ 1 Gyr. 

The spiral stellar density wave considered so far has a pitch angle oi i — 25°. 
This is typical for Sc galaxies (Kennicutt, 1981). The efficiency of outward 
radial transport of dust may depend on the openness of stellar spiral arms. 
To test this conjecture, we explored a series of models with difi'erent pitch 
angles. The dashed line in Fig. 8 shows the the dust mass outside corotation 
as a percentage of the total dust mass for the stellar spiral with a pitch angle 
i — 15°. Roughly 8% of the injected dust mass is pushed outside corotation, 
as compared to 10% for i = 25°. This indicates that the efficiency of out- 
ward radial transport of dust drops by roughly 2%. We expect that the radial 
transport of dust by spiral stellar density waves may becomes inefficient for 
very tightly wound spirals with a pitch angle i ^ 5°. Our simulations have 
also shown that the efficiency of outward transport of dust is weakly sensitive 
to the angular velocity Ogp of stellar spirals. For instance, a stellar spiral with 
flsp = 21 km kpc~-^ can transport roughly 1% less dust (see the dotted 
line in Fig. 8) than a spiral with flgp = 24 km s~^ kpc~^. We note that the 
corotation radius for Ogp = 21 km s~^ kpc~^ is at 8.8 kpc. 

The solid line in Fig. 8b shows the gas mass inside corotation as a percentage of 
the total initial gas mass. ^ It is evident that the gas (contrary to the dust) does 



^ Due to the imposed outflow outer boundary condition, we cannot reUably cal- 
culate the gas mass outside corotation because its value can be affected by the 
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Fig. 9. The ratio of the frictional timescale Tc to the dynamical timescale Tdyn as a 
function of radius at t = 1340 Myr. The soHd and dashed Hues show Tc/rdyn for the 
dynamical timescale defined as Tdyn = ^d"^ and Tdyn = |^^d^ ~ ^sp^l) respectively. 
The detailed explanation is given in the text. 

not exhibit a noticeable radial migration. A very weak inward radial transport 
(~ 1% by mass) is seen for = 24 km s~^ kpc~^, and it becomes negligible 
for Slgp = 21 km kpc~^. We attribute this difference between the gas and 
dust radial transport rates to the restoring force of the pressure gradients 
which are always present in the warm (Tg ~ 10^ K) gas. The warm gas resists 
radial redistribution by creating additional pressure gradients which act to 
restore the initial centrifugally balanced configuration. The dust, however, is 
cold and it can easily flow in the radial direction provided that the coupling 
between dust and gas is not very strong. 



To determine the main cause of the radial transport of dust, we compute 
the quantity Tc/rdyn (averaged over 0.5 kpc radial annuli). This quantity as 
a function of radius is shown by the solid line in Fig. 9 at the time when 
the gas spiral is fully developed. The coUisional timescale Tc of energy and 
momentum exchange between the dust and gas particles is defined in Eq. (4). 
The local dynamical timescale is defined as Tdyn = f^J^ — '^I'^a,4>i where fd,<^ is 
the azimuthal component of dust velocity v^. If Tc <^ Tdyn the dust is expected 
to be strongly frictionally coupled to the gas and to essentially track the gas 



gas that leaves the computational area. The dust, however, never reaches the outer 
boundary. 
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motion. Figure 9 shows that Tc is a factor of 0.01 — 0.1 smaller than Tdyn, which 
implies a strong frictional coupling betwen dust and gas. However, judging 
from the bottom panels of Fig. 7, it would appear that the dust does not follow 
exactly the gas motion in the region near corotation. The reason is that the 
dust particles have natural resonant frequencies near corotation (remember 
that corotation determines the position of resonance in gas/dust orbits). If 
the gravitational field generated by spiral structure perturbs a dust particle's 
orbit at or near its resonant frequency, then the responce of the orbit will be 
large, even when the perturbing field is weak and/or the frictional coupling 
is strong. To investigate this effect, we re-define the local dynamical time as 
Tdyn = l^d^ ^^^^1 and plot the corresponding quantity Tc/rdyn in Fig. 9 by the 
dashed line. It is now evident that the frictional coupling between dust and 
gas is much weaker near corotation. This is consistent with the trajectories 
of dust particles found in Sect. 4 - only the dust particles that are born near 
corotation show substantial radial migration. 



Another source of uncertainty in our simulations is the friction coefficient A, 
which we have adopted in the form of Eq. (4) to simplify the computations. 

As is pointed out in Sect. 2.1, the total stellar+gas surface density Eg + Eg 
(instead of simply gas surface density. Eg) should appear in the definition of 
the gas vertical scale height 



"° = 7rG(E;+Eg)' 



which may considerably increase the value of A. The uncertainity in the value 
of Eg (we do not explicitly follow the evolution of the stellar disk) works against 
this re-definition of Zq in our simulations. In order to estimate the effect that 
the increased surface density of gravitating matter in Eq. (14) may have on the 
radial transport of dust particles, we notice that in the region of interest near 
corotation the ratio of stellar-to-gas surface densities is usually in the 1.5 — 15 
range. Therefore we performed a test run in which A was defined by Eq. (4), 
but was later multiplied by a factor of 10. The efficiency of outward radial 
transport of dust increased by less than 1%. In another test, A is decreased by 
a factor of 10, and the efficiency of outward radial transport of dust is dropped 
by approximately 1%. We conclude that an order of magnitude variation in A 
has a minor effect on the efficiency of radial transport of dust. This reinforces 
a previous conclusion: the resonant region near corotation, where fl^ ~ ^^sp, 
and the coUisional drag force is dynamically less important than the spiral 
gravitational field, determines the radial migration of dust. 
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6 Dust destruction 



In this section we study the effect of dust destruction by SNe on the efficiency 
of radial transport of dust in spiral galaxies. The dust grain destruction by 
SNe is the most important mechanism for cycling the dust back to the gas 
phase. The dust destruction rate can be determined as (Dwek, 1998) 

Dd{r, (j), t) = mdest{r, (f>, t) v{r, 0, t). (15) 



Here, v{r, (f), t) is the supernova rate per unit area and TOdest(^, <P, t) is the mass 
of dust that is destroyed by a single SN expanding at location (r, 0) and time t. 
Only SN remnants with velocities exceeding 100 km s^^ can efficiently destroy 
dust grains. Therefore, the value of mdest{i^, 4>, t) is usually estimated from the 
relation nicest/ = eM^^^/Mg, where M^^^ = 3.8 x 10^ (Lisenfeld and 
Ferrara, 1998) is the mass of a remnant accelerated to velocities > 100 km s^^ 
by the blast. This relation simply states that the fraction of destroyed dust 
(by mass) within a computational cell is directly proportional to the fraction 
of accelerated remnant mass. The efficiency of dust destruction is e ~ 0.1 
(McKee, 1989). 

The SN rate per unit area is defined as 

/ T^sFR{r,(j),t - T{m)) m~" dm 

iyir,M = '^ , (16) 

/ m}~°'dm 

miow 



where Esfr is the star formation rate per unit area, a is the slope of initial 
mass function (IMF), r(m) is the lifetime of a star of mass to, toiow and TOup 
are the lower/upper cutoff masses, respectively. In practice, we approximate 
the integral in Eq. (16) by the sum with a mass discritization dm — 1 Mq. 

Observations of both normal and starburst disk galaxies suggest that on the 
scales of a few kiloparsecs star formation may be represented by a Schmidt 
law (Kennicutt, 1998) 



-•SFR 



— ^ U 2.5 X 10-^ E,^-^ PI . (17) 
yrkpcV ^ VpcV 



We modify the star formation law (17) by assuming that star formation is 
suppressed when Eg < TOm(Ej?°°^^ ^^r"")- 

Often, the instantaneous recycling approximation is used to solve Eq. (15), 
which neglects the time delay that enters Eq. (15) through the SN rate u. 
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However, we follow a more accurate approach and form stellar clusters each 
Myr according to the star formation law (17). Stellar clusters are assigned po- 
sitions and velocities drawn from the parent gas and are evolved as coUisionless 
particles in the combined gravitational potential of the halo and stellar disk. 
The local self-gravity among stellar clusters is neglected. Since each stellar 
cluster carries information on the star formation rate at the time of its birth, 
the supernova rate z/(r, and the dust destruction rate Dd{r,(l),t) in each 
computational cell (r, (p) at a time t can be easily obtained. This approach 
was successfully applied by Vorobyov (2003) to model the Ha luminosity in 
ring galaxies. We use the Salpeter IMF with a — 2.35 and lower/upper cutoff 
masses toiow = 0.5 Mq and mup = 40 Mq, respectively. 

We modify the hydrodynamic equations (l)-(3) of our gas+dust system to take 
into account star formation and dust destruction. The continuity equations for 
the gas and dust become 

BT, 

-^ + V-{E,v,)=S,-D„ (18) 



dt 



+ V-(Edt;d) = -^d. (19) 



The resulting momentum equations are 



VPg + Ed/ + 

(20) 

Ed/ - v^D^. (21) 

We have rewritten equations of motion in the form of momentum equations for 
E V which are appropriate for the system with sources and sinks. Star forma- 
tion depletes the gas reservoir of our model galaxy and is taken into account 
by the sink term Dg — Esfr in Eqs. (18) and (20). The gas is returned to the 
system with supernova explosions and quiet mass loss of intermediate and low 
mass stars. Since we do not use a multiphase description of the interstellar 
medium, the ejected hot gas is directly transformed into a warm (T ~ 10^ K) 
phase. The rate of gas ejection per unit area is determined by 



— EgVg + (vg • V)Eg'Ug = -EgV*si,s2,h - 
d 

— EdVd + {vd ■ V)Ed'yd = -EdV*si,s2,h - 



ruup 

J ^(m)EsFR(r, 0,t-T(m))m"°dm 

^g = — ^p , (22) 

/ m}''°'dm 
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where g{m) denotes the gas mass ejected by a star of mass m (Koppen and Ari- 
moto, 1991). We found it computationally prohibitive to compute accurately 
(i.e. by taking into account the time delay T(m) in Eq. [22]) the gas ejection 
rate by stellar clusters on a time scale of interest (a few Gyr). Therefore, for 
the stellar clusters older than 100 Myr, the instantaneous recycling approxima- 
tion is assumed by neglecting the time delay in Eq. (22). These stellar clusters 
instantaneously release all gas that can be produced by the 0.5 — 5 Mq stars. 
Further dynamical evolution of such clusters is not computed. 

We study the effect of dust destruction on the radial transport of dust using 
a prototype model in which the spiral density wave has a pitch angle i = 25° 
and angular velocity = 24 km s~^ kpc~^. This model was studied in detail 
in Sect. 5 where we neglected the effect of dust destruction. It was shown there 
that about 10% of dust by mass injected instantaneously inside the corota- 
tion circle can be transported outside this region by the gravitational drag 
associated with spiral density waves. We reproduce this result in Fig. 10a by 
the dashed line, and the solid line gives the dust mass outside corotation as a 
percentage of the initial total dust mass when dust destruction is taken into 
account. As is seen, supernova explosions have a minor effect on the radial 
transport of dust, decreasing the relative amount of dust outside corotation 
by only 2% when compared to the model without dust destruction. It is worth 
noting that inside corotation the destruction of dust by supernova explosions 
is efficient. Indeed, the dashed-dotted and dotted lines in Fig. 10b show the 
dust mass inside corotation as a percentage of the initial total dust mass when 
dust destruction is taken and not taken into account, respectively. Approxi- 
mately 15% of the instantaneously injected dust mass is left inside corotation 
after 2 Gyr. Given that another 10% has been transported outside corotation, 
SNe have destroyed about 75% of the injected dust mass during 2 Gyr. In 
Sect. 5 we argued that the ultimate fate of dust grains would depend on their 
place of birth. The simulations discussed in this section reinforce our earlier 
conclusion. In our model dust grains are formed inside the corotation circle at 
Tcr ~ 8 kpc. The dust grains that are born in the inner 6.0 kpc circle (by im- 
plication they constitute the majority of the total dust mass) move on elliptic 
orbits periodically crossing the spiral arms and will be inevitably destroyed 
by supernova explosions after several revolution periods. However, the dust 
grains that are born in the 6 kpc <^ r < 8 kpc annulus and sufficiently away 
from spiral arms (by implication they constitute the minority of the total dust 
mass) manage to survive the destructive effect of supernova explosions during 
many revolution periods of our model galaxy. These dust grains may move on 
orbits similar to those shown in Fig. 5 by the solid and dotted lines, circulating 
between the spiral arms (but never crossing them) and periodically traveling 
in and out of the corotation circle. We conclude that a (small) portion of dust 
grains may survive the destructive effect of SNe during many (> 10) revolution 
periods of a spiral galaxy. 
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Fig. 10. a) Dust mass outside and b) inside corotation as a percentage of the total 
initial dust mass. The solid and dashed lines show the relative dust mass outside 
corotation in the models with and without dust destruction, respectively. The dotted 
and dotted-dashed lines give the relative dust mass inside corotation in the models 
with and without dust destruction, respectively. 
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7 Summciry 



In this paper we study the outward radial transport of interstellar dust grains 
which is connected with the perturbation of circular motions caused by spiral 
stellar density waves in disk galaxies. We consider a model spiral galaxy in 
which star formation (and, by implication, dust formation) is localized within 
the corotation circle. We find that the combined action of the gravity force 
of spiral stellar density waves and the drag force between the gas and dust 
components imparts to dust particles a radial velocity component, and thus 
can evacuate dust outward. We conclude that: 

1. Dust grains that are formed inside corotation can be transported to a dis- 
tance that exceeds the corotation radius by roughly 25%. In particular, if coro- 
tation is located at rcr = 8.0 kpc, dust grains formed inside the 6 kpc ^ r < 
8.0 kpc annulus can travel to distances as large as 10 kpc in roughly 1.0 Gyr. 
This is approximately an order of magnitude faster that can be provided by 
the interstellar turbulence. 

2. A fraction of dust grains can be trapped into the interarm region, and thus 
can escape a hostile environment from supernova explosions acting mostly in 
the spiral arms. These grains form diffuse spiral arms which extend 4-5 kpc 
from the sites of active star formation. Dust particles that cross the spiral 
arms during their drift outward can survive if they do not meet strong shocks 
from SNe while passing through the arm. 

3. In the model that considers instantaneous dust injection, about 10% of dust 
by mass is evacuated outside the corotation circle in roughly 1.0 Gyr. One can 

assume that continuous injection will result in the same percentage of evacu- 
ated dust. The efficiency of outward radial dust transport is weakly sensitive 
to the rotation speed of stellar spirals, but it may be considerably reduced for 
very tightly wound spirals with a pitch angle i ^ 5°. Dust destruction by su- 
pernova explosions has only a minor effect on the efficiency of outward radial 
transport of dust by spiral density waves, which decreases the relative amount 
of transported dust mass by ~ 2%. A (small) portion of dust grains may sur- 
vive the destructive effect of SNe during many (> 10) revolution periods of a 
spiral galaxy. 

4. Our modeling shows that the dust-to-gas ratio in the region near corotation 
(±2 kpc) can vary substantially over the scales comparable to the thickness 
of spiral arms, ~ 1 kpc. Possible observations of such variations would be 
suggestive of the proposed mechanism of radial dust transport. 

In our model we neglected the fact that charged dust grains are strongly 
coupled to the magnetic field - the gyration frequency for a typical grain 
with radius a ~ 0.1//m and grain charge Z ~ 100 is of ub ~ 10~^° Hz. We 
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expect that in more realistic models, where the magnetic field is permanently 
regenerated by a turbulent dynamo, dust particles will be involved in both 
the turbulent diffusive motions and regular flow, and the overall picture will 
remain qualitatively similar to that described above. We leave this issue for 
further study. 

We should mention that the interaction of gas and solid particles through the 
drag force and the associated radial flow of solid particles have been already 
extensively studied in the context of migration of planetisimals in protoplan- 
etary disks. For example, Weidenschilling (1977) has noticed that in proto- 
planetary disks the planetesimals should migrate in the direction of increasing 
pressure (i.e. inward, for an axisymmetric disk with a declining radial density 
proflle). This is due to the fact that the gas in protoplanetary disks generally 
orbits the central protostar at the sub-Keplerian velocity due to a substantial 
contribution of pressure gradients in the gas rotation curve. The cold planetes- 
imals (which by implication have no pressure) are then forced to flow inward 
if the frictional orbital coupling between the gas and planetesimals is not too 
strong. This mechanism however may not work on galactic scales because 
the contribution of the pressure gradients to the gas rotation velocity may 
be negligible (see Sect. 2.2). In fact, the pressure gradients may even assist 
the gravity in some disk galaxies with a ring-like radial distribution of gas. 
Another interesting aspect of gas-planetesimal interaction in self-gravitating, 
protoplanetary disks was investigated by Rice et al. (2004). They considered 
disks with a flocculent spiral structure and demonstrated that the drag force 
rather than gravitational field of spiral arms plays a dominant role in the dy- 
namics of planetesimals of intermediate size, forcing them to concentrate in 
the gas spiral arms. The non-axisymmetric gravitational field of grand-design 
spirals (considered in the present paper) is by implication much stronger than 
that of flocculent spirals. As a consequence, the trajectories of dust particles 
can be considerably influenced by the spiral gravitatioanl held (Fig. 5). We 
believe that the resonant action of spiral gravitational field on dust particles 
(rather than the drag force between the gas and dust) is the main driving 
force for the radial outward transport of dust in grand-design spiral galaxies. 
This conjecture is justified by the fact that the radial flow of dust takes place 
near the corotation resonance where a weaker frictional coupling between the 
dust and gas occurs (see Fig. 9). 
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